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Abstract 

We prove that the rank 3 analogue of the tensor model defined in [arXiv:ll 11.4997 
[hep-th]] is renormalizable at all orders of perturbation. The proof is given in the 
momentum space. The one-loop 7- and /3-functions of the model are also determined. 
We find that the model with a unique coupling constant for all interactions and a 
unique wave function renormalization is asymptotically free in the UV. 
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1 Introduction 



Approaches to one of the most important problems in physics, namely the quantum gravity 
(QG) conundrum, have evolved a lot since the last two decades. The most known contender 
having some undeniable results is certainly String Theory pQ whereas alternative approaches^] 
like Asymptotic Safety scenario [2], Noncommutative Geometry [3J, Dynamical Triangula- 
tions [4J and Loop Quantum Gravity (LQG) [5] have also drawn a lot of attention of theorists 
from the mid 90's (see for instance [6J [7J [8J |5j) when, in fact, the String's revolution enchanted 
most of the physics community. 

In the meantime, another framework starts to build up around the Sakharov's idea (1965) 
of an "emergent" theory of gravity (see [9] for a review). Mainly "emergent" refers to a 
phenomenon which is only induced and not fundamental. The analogy with hydrodynamics 
as emerging from laws of molecular physics is commonly used as an illustration. This idea, 
somehow rooted in condensed matter physics and statistical physics, suggests only, but very 
originally, that the quantization of the spacetime background metric should be addressed 
in some independent way than the ordinary study of fluctuations around the flat metric 
configuration. In fact, it is well known that the latter type of quantization leads inexorably 
to a non renormalizable quantum field theory using Einstein-Hilbert action [TQ1III]- So far, 
results on the QG domain are various and are greeted with more or less success. Nevertheless, 
the idea of an emergent spacetime was welcome as pertinent and still perpetuates through 
the years [T2] . 

Coming back to the mid 80's, the so-called random matrix models, initially motivated 
by string theoretical considerations, become on their own a concrete proposal for statistical 
models for QG in 2D [13]. Soon after, they were followed by their tensor analogue for 
D > 2 QG [14|-[16j. However, the latter works experienced difficulties because one of the 
main analytical tool, namely the 1/N expansion, was crucially missing. In last resort, only 
numerical results can be properly achieved. All these endeavors appropriately realize the 
idea of an emergent theory of gravity (in short, in 2D, random matrix models are theories of 
random triangulations of surfaces and summing on such triangulations amounts to sum over 
geometries of these surfaces. These models possess a phase transition towards a "continuous 
phase" geometry the so-called conformal geometry coupled to Liouville gravity [13J). 

The framework of tensor fields for addressing the QG issue was correctly stated in [T7] 
for a D-dimensional lattice version of 3D Euclidean gravity. The prime idea was not so much 
focused on the understanding of an extended version of statistical analysis of matrix models 
for tensors but rather to introduce a discrete version of a simplified and "flat" (BF) model 
of gravity valid in dimension higher than 2. Substituting the integral over all D-dimensional 
manifolds by a sum over all D-dimensional simplicial complexes, the ensuing models generate 
perturbatively all D simplicial complexes providing the tensor analogues of the successful 
matrix models. Fields here are defined over several copies of an abstract Lie group which, 
by Fourier transform, yield tensor components. Depending on the dimension, to each tensor 
is associated a simplex and vertices provide fusion rules and exchange of momenta as in 
ordinary quantum field theory. A key point of this formalism is that these tensors should 
be constrained to satisfy specific rules of invariance enforcing the flatness condition of the 

include Supergravity and M-Theory within the String approach as theories invoking extra symmetries 
or dimensions. 
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gluing of the simplices at the vertices. A generic link was later found between these lattice 
models and another fast-growing field: the spin foam models which embody the covariant 
version of the LQG program [TBI [19] . Thus, was born a new line of investigation for QG, 
the so-called group field theory (GFT) [201 l2Tj . The partition function of GFT sums over 
both topologies (dictated by the topology of Feynman graphs as simplicial manifolds) and 
geometries (encoded in the amplitude of these graphs) of a given manifold. GFT claims to 
be a theory of quantization "of" spacetime itself [21]. This leads, once again, to the problem 
of obtaining an emergent spacetime at some proper continuum limit. 

The GFT framework appears very appealing for performing quantum field theory compu- 
tations and, naturally, the question of its renormalizability was first systematically addressed 
[221 123]- Many other interrogations arose concerning the GFT formalism. "How to control 
divergences ?" and "why the most important contributions in the partition function would 
be of the form of a "large" and "smooth" spacetime like the one, we experience ?" were 
frequently asked. In fact, all these interrogations could have been only satiated by finding a 
1/N expansion for tensor models. 

The renormalization program for GFT made its first steps and notable facts have been 
sorted out concerning power-counting theorems of both ordinary GFTs and more involved 
GFT models in relationship with 4D QG [24j-[28j. In the meantime, a drastic improvement 
concerning the topology of simplices dual to the Feynman graphs has been proved by Gurau 
by introducing the colored version of these GFTs [29] [30]. In |3~T] . it appears clear that 
relevant operators of the Laplacian form [32] were missing in the 3D Boulatov GFT action 
[T7] and, so, should be added in all GFT actions before discussing of their renormalizability. 
The renormalization analysis for these type of GFTs proves to be a computational challenge. 
Concerning the symmetry analysis of GFTs, some efforts have been made in order to show 
that they possess a quantum group symmetry [33J, that they satisfy peculiar Ward identities 
[3l] and that, associated with translation invariance, there exists a conserved classical energy 
momentum tensor for the colored GFTs [33] . 

The real breakthrough of the story occurs with the major discovery by Gurau of the 
analogue of 1/N expansion both for GFTs and independent identically distributed tensor 
models (see [36J [38J and more other references therein, |39j for a short review of the subject 
and [40J for a complete overview of the subject and more developments based on this 1/N 
expansion). Dominant graphs were identified [UJ and are associated with simplices of the 
sphere topology. Important developments followed: the generalization of the Witt-Virasoro 
algebra (without central charge) for infinite dimensional tree algebras [42] . the tensor gen- 
eralization of the Ising model in any dimension [13], the determination of the universal 
character of random tensor models generalizing the Wigner-Dyson law for random matrices 
[Hj. Furthermore, based on the above developments and for the first time, a 4D QG model, 
even though simplified, has been found renormalizable at all order of perturbation theory 
[45] . Let us readily mention that it is not clear in which sense General Relativity could be 
recovered from the said model (this is why one calls it simplified). However, it is not hopeless 
to encounter coordinate invariance and more geometrical properties of the background space 
emerging through the rich structure of that model [16] . 

The idea of defining a renormalizable and emergent theory of gravity using the tensor 
approach has matured through years [16] and belongs to a particular pattern extending pro- 
gressively the renormalization group (RG) analysis from local graphs like in the <p\ model, to 
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the vector cases for condensed matter systems, then to the matrix cases like in noncommuta- 
tive field theory [17]- [51] incorporating already nonlocality and, now, to higher dimensional 
tensor theories. The model of [15] is nonlocal and defined over copies of the compact group 
£7(1) hence does not have infrared (IR) divergences. It is built by first integrating four colors 
(out of five) in the colored model [29J and choosing the last color to be dynamical with re- 
spect to [31J. Considering the 1/N expansion, it is known that only certain contributions in 
the effective action called "melonic" (this terminology refers to [UJ) will be not suppressed 
in power-counting. It remains to truncate the effective action for obtaining relevant and 
marginal operators of the melonic kind. After analyzing the divergence degree of an arbi- 
trary graph at a given order of perturbation theory, only graphs determined by particular 
boundary data but having also a regular internal structure turn out to be divergent. This is 
the "generalized locality" principle for tensor graphs and only those graphs should be renor- 
malized by standard interpolation moves on external data. Interestingly, a peculiar anomaly 
arises in the expansion. The authors of [15] interpreted this term as matter defects appear- 
ing in that 4D gravity model. Renormalization group (RG) flows have been not computed 
for the model and these deserve indeed to be understood. As it can be easily realized, the 
combinatorics of such a model becoming very involved should be cautiously handled. The 
present contribution is a step towards the computation of the RG flow of that 4D QG. 

In this paper, we address the RG flows of coupling constants of a tensor model in 3D 
related to the above 4D model. This model has been briefly outlined in [15] but its renor- 
malizability has been not yet proved. This is an initial property that we need to investigate 
before computing any flow. We prove that this 3D tensor model over U(l) 3 is renormalizable 
at all orders of perturbation. The proof is thoroughly done in the momentum space. This is 
in contrast with the renormalizability proof of the anterior 4D model which was performed 
in the direct space. In the momentum basis, the interaction remain of the same form of as in 
the direct space yielding factorized graph amplitudes somehow more suitable to perform the 
different optimizations occurring in the multiscale analysis. Hence, using a different basis 
also convey to a new perspective on these tensor models in general. In a second part of this 
paper, we compute the one-loop 7- (governing the RG flow of the wave function renormaliza- 
tion) and /3-functions (governing the RG flow of interaction coupling constants) and analyze 
the RG flows of the different coupling constants. We emphasize that it is not necessary to go 
beyond one-loop computations since the leading order corrections determine the RG flows 
of coupling constants if the /3-function is not vanishing at this order. The model obtained 
by merging all coupling constants to a unique one which is somehow the most natural model 
proves to be asymptotically free in the UV. This feature might be related to the universal 
Gaussian behavior of these tensor models discussed in [H]. 

The plan of this paper is as follows: the next section defines the model and states our 
two main results. Section 3 is devoted to the multiscale analysis and the achievement of 
a crude power-counting theorem which will be dissected in Section 4. The renormalization 
of primitively divergent graphs will be performed in Section 5. The calculations of 7- and 
/3-functions and RG flows of coupling constants are detailed in Section 6. Section 7 gives a 
summary of our results and an outlook of this work. An appendix gathers further details on 
some results used through the text. 
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2 Rank 3 tensor model over U(l) 

Let us consider four complex three rank tensor fields over the group U(l), f 8, : U(l) 3 — > C. 
The index a = 0, 1, 2, 3 is called color. In Fourier modes, (p a can be expanded as 

^(9i,92,9s) = £ ^./^e^e^, t e [0,2tt) and [ Pj ] = ( Pl ,p 2 ,p 3 ), (1) 

where the group elements g% = e tdk G (7(1) = S 1 . From now on, we write (fl 23 := ( P^ P2P3 
and define a theory in the momentum space with a first kinetic part regarding the colors 
a = 1, 2 and 3, 

^ iW ^t&^ (2) 

p-j a=l 

where the sum in pj is performed over all momenta values, for j = 1, 2, 3. We do not assume 
any symmetry under permutation of arguments of these tensors. 

The ordinary colored theory [29] is defined by an interaction which can be read off in 
momentum space as 

S mt = A £ Vl23 ^345 ^526 ^641 + * Yl ^123 ^345 ^526 ^641 > ( 3 ) 
Pj Pj 

the parameters A and A being the coupling constants. 

We emphasized here an important point concerning these tensor models. The momentum 
space is in "exact duality" with the direct space in the following sense: up to an inessen- 
tial constant (a power of 2tt coming from spatial integrations), interactions in direct and 
momentum spaces share exactly the same form. This is contrast with other local theories 
or even nonlocal theory like in noncommutative field theory [JT] for which the vertex pos- 
sesses a delta function of momentum conservation. Although nothing prevents to perform 
the renormalization in the direct space, the main point for switching in momentum space is 
that graph amplitudes get factorized in a different and useful way as we will see. 

The polar point here is to take a particular kinetic term with respect to the color such 
that 

3 

S kin '° = £ ^? 23 ( £ a. \ Ps I + m) ^? 32 , (4) 
Pj 3=1 

where \p s \ denotes the absolute value of p s and a s are positive wave-function coupling con- 
stants. Hence the field of color is assumed to be propagating]^] Note that also, the model 
described by Q is slightly different from the direct 3D analogue of [13] . The latter could be 
only reproduced by fixing all a s to a given value. 

We use the same procedure of [12] which mainly performs an integration of the partition 
function with respect to all colors save one. One gets an effective action for the last tensor 

2 Remark that a direct space formulation corresponding to such a momentum space kinetic term can 
be defined by a reduced operator acting on each strand as —idgcj){d) = X) P ez |p|<^(p)e zpe , where <j)(p) is the 
Fourier mode of cf>. 
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Figure 1: Propagator. 



(p° in the form 

Z= [d»c[v ]e- slat '°, (5) 



where dfic[f°] stands for the Gaussian measure with covariance C = (J2 s a s\Ps\ + m ) 1 
(represented in Figjl]) and the effective interaction finds the form 

&* fi = Y #^iV /(p ' D) -^ w(e) Tr B [^V ]; (6) 
^ Sym(i3) 

the sum in B is performed on all "bubbles" which are connected vacuum graphs with colors 
1 up to D = 3 and p vertices; f(p,D) is a function of the number of vertices and the 
dimension; oj{B) := J2j9J i s the sum of genera of sub-ribbon graphs called jackets J of 
the bubble, and Trgf^ ^ ] are called tensor network operators. Graphs with u(B) = are 
called melons [ID] and non melonic contributions defined by u(B) > become suppressed 
from ([6]). Thus, it is sufficient to only focus on the melonic sector of the theory. As in [42], 
one attributes a different coupling constant to different tensor network operators and simply 
rewrites (omitting the color index 0) 

We truncate the above series and will consider only relevant to marginal terms guided by 
renormalization conditions. Thus, in the following, we will consider the effective interaction 
terms made of monomials with order at most four: 

L Ppi,P2,P3 ( Pp 1 ',P2,P3 L Pp 1 r,p 2 ',P3' ^pi^'iPs' + P erm utations. (8) 

Pj 

The last term "permutations" means that we include also other terms by performing a 
permutation over the six momentum arguments (see Fig. [2]). 

By introducing the UV cutoff A on the propagator which becomes C A , the action and 
partition function of our model are then defined by 

" d^]e~ s \ (9) 



S A — \ A 1S4 + CT^ S2-1 + CT A 2 - 2 Ss;2;2, Z — I • 

s=l J 



where 

s 2]1 = J2n Pj mpjh s *w = '52<p\pj]{\p*\) ( Pipj]> ( 10 ) 



pj pj 

where the symbols CT's are coupling constant counterterms defined by the difference between 
the renormalized and bare couplings. More precisely, CT^ is a mass counterterm and CT S A 2 . 2 
is a wave function counterterm. 

The main results of this paper are given by the following statements: 



Figure 2: Vertices of the type V4. 



Theorem 1. The model defined by ([9]) is renormalizable at all orders of perturbation theory. 

Theorem 2. The model obtained from ^ by identifying a s = a is asymptotically free in the 
UV direction. 

3 Multiscale analysis 

We begin with the multiscale analysis which will lead to a prime (or crude) power-counting 
theorem. From that theorem, we will be able to identify the "dangerous" graphs for which 
the renormalization program should be performed. The first step is to find the behavior of 
the propagator with respect to high and low momentum scales, then, using this result, we 
should find an optimal way to bound graph amplitudes in the most general manner. 



3.1 Propagator bound 

Consider the kinetic term of the action. The kernel of the propagator is 



3 _ 3 



J Ps,P' s 

S=l S=l 



the notation C(\p s ], \p' s ]) referring to C(\pi,p2,Pa], [p'^P^P's})- 

Introducing a Schwinger parameter, we get the integral form of the propagator as 

POO 

C(\p,]M)= dae- a[ ^ aM+m] . (12) 
Jo 

In the following developments, only is needed the small distance or high momenta behavior 
of that propagator. We do not have the infrared divergence problem from the compactness 
of U(l). This is reflected in the direct space (Fourier transform of the above), for instance, 
by the fact that the mass can be put to m = and this will cause no difficulty with the 
zero mode of the propagator. Working in the momentum space, in full generality, it is 



better consider a non vanishing mass otherwise low momenta will cause a divergence in (12). 
Nevertheless, as mentioned earlier, this situation is actually not relevant for the remaining 
analysis which solely focuses on the UV sector of the theory. 
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The next stage is to introduce the slice decomposition of the propagator in the form: 



C = ^TCi i = 1, 2 • - • (13) 

8=0 

poo pM~ i+1 

C = dae~ a[ ^ asM+m] , C t = dae- a[ ^ asM+m \ Mel (14) 

Jl J M~* 



The following proposition is direct 

Lemma 1. For all i e N, there exists a large constant K > such that 

Ct < KM- i e- M ~ i ^=i a ° M+m l (15) 
The ultraviolet cutoff can be imposed by summing the slice index % up to a large integer 



called A in the slice decomposition (13). Thus 

A 



c a = J2 c ^ ( 16 ) 



i=0 



and the ultraviolet limit consists in taking A — » oo. One calls Cq and C\ the IR and UV 
propagator slice, respectively. For simplicity in the following, we omit the superscript A. 

3.2 Optimal amplitude bound: Prime power- counting 

Let Q be a connected amputated graph with set of vertices V (of any kind for the moment) 
of cardinal V and £ set of lines of cardinal L = \C\. The bare amplitude associated with Q 
is of the form 

n 

Pv,s t&C veV; s 

where /^[A, CT 2j i, CT S . 2; 2] is a function of some product of coupling parameters, S(Q) is 
a symmetry factor, p v (t) }S are momenta involved in the propagator which should possess a 
vertex label v(£) hooked to a given line I and a strand label s; p VtS are the same momenta, but 
now, involved in the vertex which should bear both vertex v and strand s indices; 5 Pv stPv , 
is the Kronecker symbol (we keep here formal notations but, while dealing with a category 
of vertex V4, V2, this symbol might depend on the vertex type and the number of strands s 
could vary from one category of vertices to the other; however, as we will see below, after 
some momentum summations, we do not need to track all these indices to get a crude power 
counting); \i = i 2 , • • • ,iq) is called momentum assignment and gives to each propagator 
of each internal line i a scale i# G [0, A]; the sum over /1 is performed on all assignments. 
From the point of view of the effective series expansion [52], the function / M collects all 
the effective couplings corresponding to the attribution fi. Given an amputated graph, we 
simply have external vertices (where test functions or external fields can be hooked). By 
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convention, we fix all external line scales at i ext = — 1. Ag.^ will be the core quantity and the 
sum Ag = ^2J_f l j,(^,CT___)/S(Q)]Ag. l j j can be done only after renormalization in the standard 
way of [52] . 

We would like to perform the momentum sums in the p V)S in an "optimal" way. For that 
purpose, let us quickly review what should be expected from the ordinary theory. Given \i and 
a scale i, we consider the complete list of the connected components G*, k = 1, 2, . . . , k(i), 
of the subgraph Qi made of all lines in Q with the scale attribution j > i in \x (note also that 
the meaning of i and k are radically different even though we keep simple notations for these 
quantities when writing Gf). Such subgraphs are called high or quasi- local and are the key 
objects in the multiscale expansion [52]. There is a partial (inclusion) order on the set of G\ 
and Qq = Q . The abstract tree made of nodes as the G\ associated to that partial order is 
called the Gallavotti-Nicolo tree [S3]- Q is the root of that tree. To an arbitrary subgraph 
g, one assigns two quantities: 

i g (fi) = iniii(fi) , e g (/j) = sup ii(^). (18) 

'^fl l external line of g 

A subgraph g can be viewed as G% for a given /i if and only if i g (fi) > i > e g (/z). In the 
direct space formalism, one considers a spanning tred^T of the graph Q. Associated with T 
are position variables that one integrates to give decay factors to the product of propagator 
lines. The key point is to optimize the bound over spatial integrations by choosing the tree T 
to be compatible with the Gallavotti-Nicolo tree. This is achieved by taking the restriction 
of T to any G\ such that Tf is again a spanning tree for G\. 

Coming back to our situation and for simplicity, we assume that no wave function coun- 
terterm appears in the graph. Adding them at the end will be an easy task. One notices 
that the vertex operator is a product of delta functions and hence Ag factorizes in term of 
closed or open strand line that we call "faces". Let T be the set of faces of Q. It decomposes 
in a set J 7 -^ of closed (or internal) faces and another set J-" e xt of faces connected to external 
vertices. We have 



and 



\F\ — F — F- mt + F ext , F int — iJ^nt |, F ext — |.F ext |, (19) 

i^i < ^ m n M ~"En[n^ s <J e " M " [Esasksi+mi 

£&C q s fe£ s=l 

< K' n Y[M- i ^HIl e ~ M ~ Hafl9f ^ (2°) 

tec q } fe^eef 



where is the wave function quantity now bearing a face label /. The bound (15) has 



been also used in order to get (20). The following cases may occur: 

(i) The face / G J 7 ^ , then the face amplitude is £ e -(£* 6 / M ^WM. Given i = 
mmifzfii, this amplitude can be optimized as, up to some constant S, 

^V(£ <6/ A/- ! 0a/k/l < ^ e -Af-Wkl = §m { + 0(M-*). (21) 
if i 



3 This is by definition a subgraph formed by lines passing through all vertices without forming closed 
loops. 



s 



(ii) the face / is open, then all sums in q s can be performed and yield 0(1). 



Hence, in the above amplitude (20), only terms involving J 7 ;^ should be taken into account. 
We obtain 

\Ag.J < K' n J[M-^Y, II e-Eie/^Oa/lff/l 
l£C q f fen* 

h 

< K' n YlY[M 1 Yl K"M il f, (22) 
iec i=i fefint 

where If is a strand of / such that i\. = inf; e / %%. It can be therefore inferred that 

\A g .j < K ,n K" F ^ yi n m_i n ri M 

e&c i=i /e^nt t=i 

< K' n K" F ^n ii m- 1 yi n m 

< K >n K » F ^ y[ m l g ) n n n m 

(i,fc)eN 2 (j,fc)6N 2 /eJ'int nG, fc iftfnG* 

< tf/n^z/JU JJ M -i(G*)+ J F int (Gf) ) (23) 

(i,fc)eN 2 

where L{G\) and F int (Gf ) denote the number of internal lines and internal closed faces of 
G\. In the last step, during the bound optimization, one uses the fact that If 6 /fl G# gives 
a single strand the scale of which is the minimum among the scales of the lines which occur 
in / n G\. Note that the face /nG- can be an open face. In such a case, If being the 
minimum of the scales of the closed face / cannot occur in the / D Gf, hence the product in 
If G / fl G* is empty. Thus, the last product yields nothing but the number of elements of 
the internal closed faces of G\ (when the closed face / becomes totally embedded in the G\ 
such that If G G\\ see Figj3|. 

Compounding all constant factors coming from the momenta sums, we get a bound of 
the graph amplitude at a given attribution fi as 

\Ag ;fl \ <K n Y[ M- L W +F **W\ (24) 
(M) 

where K is some constant and n is the number of vertices of the graph. We recall that 
the above amplitude is assumed to be without wave-function counterterms. Adding these 
counterterms in the form of V s . 2 vertices, the following statement is straightforward: 

Lemma 2 (Prime power-counting). For a connected graph Q (with external arguments in- 
tegrated versus fixed smooth test functions) , we have 

\Ag.J < ir n JjM^ (G ^, (25) 

(i,k) 
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fl={sl, S3, s4} f2={s2, S3', s4'} 

[ 12 Sf2AG 1 1 5 l 12 (l,15)=^ 

l E «f2fl(4 l (2,12)=jzf 

! 12 Sf2AG 1 \ l n (2,U)=^ 

l 12 €f2AG 1 1 l, 2 (l,10)=^ 

[ 12 <=f2AG ! 1 I (l,9)={s4} 

l 12 6f2AG„ 1 l o (l,0)={s4} 



Figure 3: A graph with a multiscale expansion (LI is at scale i = 15, etc.) 
and its Gallavotti-Nicolo tree; the face fl (in red) is external and the face /2 
(in green) is internal. If is the strand element at given scale i which should 
optimize the bound. 




Fi gure 4: A graph Q (left) and its color extension Q co \ox (right) in simplified 
notations: each line of color a = 0, 1,2,3 corresponds to a propagator in the 
colored theory, i.e. J d \ic (</?) <f 'm^m' anc ^ vertices are defined by Eq.Q. 

where K and n are large constants, 

3 

MGt) = ~L(Gi) + F mt {G k t ) + £ VUG*)- (26) 

s=l 

The quantity Ud(G) is called the divergence of the graph Q and provides also the power- 
counting. In the sequel, we aim at analyzing this quantity for a general connected graph 

Q. 

4 Analysis of the divergence degree 

The analysis of the divergence degree is made in two steps: the first by translating Ud{G) in 
topological quantities and the second by refining the obtained result in a convenient manner. 
From these two steps, we provide an exhaustive list of divergent graphs. 

4.1 Power- counting in topological terms 

We recall some definitions (see [41j[54j[45j): 
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J J 



Figure 5: The jacket J (0123), ribbon subgraph of the colored graph Q co \ or 




eg dg dg = j 

Figure 6: The boundary dQ of Q and its rank two or ribbon structure. 



Definition 1. Consider g a 3 dimensional graph. 

(i) The colored extension of g is the unique graph g co ior obtained after restoring in g the 
former colored theory graph (Definition 1 in [J^S], here illustrated in dimension 3 in 
FigQ. 

(ii) A jacket J of g color is a ribbon subgraph ofQ co i or defined by a cycle (Oabc) up to a cyclic 
permutation (see Definition 1 of [J^l] and an illustration given by Fig\^. There are 3 
such jackets due to the dimension 3. 

(Hi) The jacket J is the jacket graph obtained from J after "pinching" that is the procedure 
consisting in closing all external legs present in J (see Section 3.2 of [5$ , for the 
general definition of "pinching" for external strands; here the pinching of a jacket is 
illustrated in Fig^. It is always a vacuum graph. 

(iv) The boundary dg of the graph g is the closed graph defined by vertices corresponding 
to external legs and by lines corresponding to external strands of g (see Section 3.2 of 
[5$ as well as a drawing in Fig^fy. Here, it is a vacuum graph of the 2 dimensional 
colored theory, hence a ribbon graph. It is also its own and unique jacket. 

Some comments are in order. In general, the boundary graph of a tensor graph of rank D 
is a closed tensor graph or rank D — 1. The notion of jacket for the boundary graph can be 
defined as in the ordinary situation as a cycle of colors. The fact that the boundary graph 
is denned as a closed graph induces immediately that all jackets are closed. 

Let g be a graph. Let V4 be its number of ip 4 vertices, V2 the number of vertices cor- 
responding to mass-counterterms (p 2 , V' s . 2 the number of vertices corresponding to wave 
function counterterms |p s Iv^ 2 . The graph possesses also L number of lines and N eXt external 
fields. Using the above notations, the following statement holds: 

Theorem 3. The divergence degree of connected graph g is given by 

u>d(G) = ~V 2 - l -{N ext - 4) -Y^gj + gog - (C dg - 1), (27) 

j 



11 



where the sum is performed on all jackets J ofQ co i or , gj is the genus of the jacket J associated 
with J, gsg is the genus ofdQ and Cgg is the number of connected components of the boundary 
graph dQ. 

Proof. There is a relation between the numbers of lines, of external legs and of vertices for 
Q: 4V4 + 2(V 2 + Yf s =i K'2) = 2L + iV ex t . Concerning £ co ior, its number of vertices Vg color and 
number of lines Lg coloi satisfy 

3 1 
Vg coloI = m + 2(V 2 + K2I Lg coloi = L + L int] g coloi = -(AV Qcolm - 7V ext ), (28) 

s=l 

where L int; g color denotes the number of internal lines of Q co \ or which do not appear in Q. 
Fg coloi , the number of faces of £ C oior, can be partitioned in the number F of faces of the initial 
graph but also in additional faces F int .g color proper to the colored graph and not appearing 
in Q. We write 

^color = F + Fint; £; co i or - (29) 

There are 3 jackets in Q co \ or . Each face of the graph Q co \ or is shared by 2 jackets. Summing 
over the jackets, one therefore has Ylj Fj = 2Fg color . Meanwhile, concerning number of 
vertices and lines, we have Vj = Vg colOT and lines Lj = Lg coloi , respectively. 

The next stage is to pass to topological numbers associated with the graphs. The Euler 
characteristic of a ribbon graph can be only defined if the graph is closed. This is the reason 
why we need to consider pinched jackets J. After closing all external half-lines of open 
jacket graphs J, we are in presence of closed ribbon graphs J for which the above topological 
number is well defined. For the resulting jacket J, we have 

Vj = Vj, Lj = Lj, (30) 

hence J has the same number of vertices, the same number of lines as J, but a different 
number of faces. The number Fj of faces of J finds the decomposition Fj = F int .j+ F cxt .j, 
where F int .j is equal to F int; j the number of internal faces of J and F ext .j is the number of 
additional closed faces entailed by the pinching procedure. 
The Euler characteristic of J affords 

F int .,j+F e ^j = 2 - 2gj-Vj + Lj. (31) 

F int .j can be further divided into F- mt .j.g, the number closed faces belonging to Q and 
F... j. r , the number of closed faces belonging to £? C oior and not to Q: 

HIT, , J , y color 

^int ; J = F int -J-g + F- mt . J. g coior . (32) 

Summing over all jackets, we have 

Yl ( F int ; J; Q + F int ; J; © color + F ext ; j) = 2F mt ;Q + 2F int . + F ext ; J ' ( 33 ) 

J J 

where, we recall that F- mt .g coloi is the number of faces issued from <? co ior and not appearing 
in Q and F mt -g is the number of closed faces of Q. 
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The quantity F int . g color can be computed explicitly: each ip 4 vertex contains 4 internal 
faces coming from the coloring (and not present in Q) and each (p 2 type vertex contains 3 of 
them. Then, one has 

3 

*i=ti 5 coloI = 4V 4 + S(V 2 + K ;2 )- (34) 

8=1 



Besides, using (28), we have 



J2 (~Vj + Lj) = 3(Vg coloI ~ \N cxt ) = 3 UV 4 + 2(V 2 + £ 

J V 8=1 



(35) 



Summing over J in (|31|), and using the above relation, we get from (J33|) and (|34|) 

1 
2 



FuA;Q — 2V4 — J-^Vext 



3 ~ 5Z#J~ 9 5Z F ext;J- 



(36) 



We re-express ^2jF exk .j i n terms of topological numbers of the boundary graph dQ. The 
boundary dQ is defined from its number of vertices Vgg and lines Lqq such that 



Vqq — N ext , Lqq — F ext . 

The external legs of the initial graph Q have 3 strands and we have 

3iV ext = 2F cxt . 



(37) 



(38) 



The boundary graph is a closed ribbon graph, thus, dQ has a single and closed jacket, itself. 
Note that a boundary graph may have several connected components. We get from the Euler 
formula 

2Cdg - 2g d g = Vqq - Lqq + F d g, (39) 

where Fgg and Cgg are, respectively, the number of faces and of connected components of 



dQ. It is simple to deduce from (37) and (38) that Lqq — Vgg = N ext /2 and, from the latter 
and the above Euler formula, one recovers 



F dg = 2{C dg - 1) - 2g dg + 2 + ^N ext . 



We remark that 



$^ F ext ; J - F 9G, 



(40) 



(41) 



because each face of the boundary graph is uniquely represented in a unique J. Indeed, we 
recall that a face of the boundary dQ can be represented by a color triple (Oafe); a face (Oab) 
belongs to a jacket J if J is the form (Oacb). A jacket being a cycle this implies that c should 



be fixed (See FigjTj). Then evaluating ( [41] ) via (40), and inserting the result in (36), it can 
be inferred 



F in t-g = 2V4-jN mt +3-J2gj-h2(C g g-i)-2g gg + 2+}-N ext ] 
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Figure 7: A face (highlighted) of dQ (boundary graph of Q Figj4]) labeled by 
colors (013) and the unique pinched jacket J (0123) which possesses this face. 



2V 4 - AU + 2 - 9J ~ (Cag - 1) + 



9ag- 



(42) 



Inserting the latter in u d {G) = -L + F hit , g + V{, where L = (1/2) [4 V 4 + 2(V 2 + £j =1 1/', 



iVext], we finally get (27) which achieves the proof of the theorem. 
The quantity 

~^29j + 9dg- (Cog - 1) 



□ 



(43) 



should be analyzed in detail since, mainly, the classification of divergent graphs relies on its 
behavior. The next section is devoted to this study. 



4.2 Bounds on genera 

This section undertakes the study of the quantity — 9 J + 9QQ which turns out to be the 
central object capturing the behavior of the divergence degree. From that result, we will be 
able to classify the primitively divergent graphs in a next stage. 

The strategy here follows mainly the same adopted in [15]: we perform a sequence of 
0/c-dipole contraction [36j E3 EE] of a given graph and scrutinize the genus change under 
such moves. 

Let us give a flavor of the following combinatorial analysis before starting it. A fc-dipole 
contraction on a colored vacuum graph [35] generalizes to the tensor case, the so-called line 
contractions along a tree for matrix ribbon graphs. Performing such a sequence erases all 
bubbles with melonic structure and one gets another graph which can be called a Filk rosette 
[57] of the tensor kind. The degree of that graph uj{Q) = J2j9J * s the same as the degree 
of the initial graph since the genus of each jacket along this sequence of contractions is 
preserved. 

In the present case (and also in [IS]), starting from a colored theory with colored graphs, 
we do not perform arbitrary A;-dipole contractions but, at first, fc-dipole contractions involv- 
ing all colors but 0. The result is a graph with the same degree as previously claimed. In 
particular, erasing all melonic contributions yields exactly a graph of our starting theory. 
In a second step, we start what we call a Ofc-dipole contraction involving now the last color 
(a precise definition will be given below). Under such a contraction, the graph may or 
may not change of degree. We should analyze in detail the consequence of performing the 
contraction. The main points revealed in [45J were: (1) the boundary graph coincides with 
the resulting graph obtained after any full Ofc-dipole contraction in an arbitrary sequence 
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02-dipole 03-dipole 
Figure 8: Ofc-dipole contractions. 



and removal of external structure of a given graph; (2) under a single Ofc-dipole contraction, 
the genus of a jacket never increases and (3) the sum over all jackets of differences between 
genera of these jackets before and after contraction is always bounded from below by a fixed 
constant. 

In the present context, we actually expect the same properties for graphs. 

Definition 2 (0A;-dipole and contraction). We define a Ok-dipole, k = 0, 1, 2, 3, as a maximal 
subgraph of Q co ior made of k + 1 lines joining two vertices, one of which of color 0. Maximal 
means the Ok-dipole is not included in a 0(k + l)-dipole. 

The contraction of a Ok-dipole erases the k + l lines of the dipole and joins the remaining 
3 — k lines on both sides of the dipole by respecting colors (see Fig^. 

We summarize in the following proposition a basic fact about the boundary graph (the 
proof is totally similar of that of Lemma 3(Graph Contraction) in [15] ) 

Proposition 1 (Graph contraction). Performing the maximal number (4V4 — N ext )/2 of 
Ok-dipole contractions on Q co i or in any arbitrary order and erasing the external legs and the 
remaining open faces of Q leads to the boundary graph dQ . 

We restrict the analysis to the only significant situation of graphs without any two-points 
vertices of the kinds V 2 and ^'.Jj and aim at studying the change in genus of a single closed 
jacket under a dipole contraction. We write 

Scolor -> ^color, J ~+ J', 9 J ~+ 9j>- (44) 

The analysis of the sum of differences between genera before and after the contraction will 
be a corollary of that result. 

Consider a colored connected graph Q co \ OT1 a fixed 0/c-dipole and the contracted graph 
G'coIot, which may or may not be connected. After a contraction, the numbers of vertices and 
lines meet the following modifications: 

V ->■ V = V - 2, L^L' = L - 4, (45) 

whereas the number of connected components may change from c = 1 to d < 3. 

We can first analyze graphs for which one performs a full and maximal contraction of all ip 2 type vertices 
into a single line. After, the case including such vertices can be derived by re-inserting ip 2 vertex chains from 
that point. 



15 



The other ingredient entering in the definition of genus is the number of faces. The change 
in faces can be handled by another notion of pair types associated with dipole contractions 



see Appendix A.l). The key observations are summarized by the following statements: 



Lemma 3 (Decreasing genera). Given a pinched jacket J, we have 



(46) 



Proof. See Appendix A.l| 

Lemma 4 (Genus bounds). We have 

gj > gag 



3g dg e N. 



(47) 



Proof. A complete sequence of Ofc-dipole contractions on the graph Q ^yields the graph 
dQ (see Proposition [l]). Using Lemma ij any genus of any pinched jacket J decreases along 
that sequence of contractions. Finally, the pinched jacket coincides with the boundary graph0 
itself dQ, which proves (47). The factor of 3 appearing in the second statement of (47), come 



from the fact that there must exist at least three jackets J (in the sequence) generating the 
boundary graph which possess a greater genus. Indeed, if (abc) defines the boundary graph 
cycle (obtained after contracting the color 0, such that a, b, c ^ 0), then three jackets (Oabc), 
(aObc) and (abOc) contract back on (abc). Some of the "ancestors" (placed upstream with 
respect to the sequence of contractions) of these three jackets or the three jackets themselves 
should possess a genus strictly greater than g$g. 

□ 



4.3 Divergent graphs 

We always assume that dQ ^ hence Cgg > 1. Furthermore Cgg < N ext /2, because each 
connected component must have at least a non zero even number of external legs. 



After an analysis (the details of which are collected in Appendix A.2), the following table 
gives a list of primitively divergent graphs: 





v 2 


gag 


Cog - 1 


Ejgj 




4 

















2 














l 


2 


1 














2 











i 






Table 1: List of primitively divergent graphs 
We emphasize the fact that the anomalous term of the form (J <p 2 ){f <p 2 ) occurring in 



four dimensions [45J does not appear here. One could claims that in higher dimensions, more 
anomalous might be present. From the above table characterizing all dangerous contribu- 
tions, we are in position to address the renormalization program for these amplitudes and 
this is the purpose of the next section. 



5 This results from Proposition [T] since after a full sequence of contractions Q — > dQ, any pinched jacket 
of Q should have a projected in dQ. However dQ is unique so that the conclusion is immediate. 
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5 Renormalization 



The renormalization procedure of primitively divergent graphs could be implement in the 
way of jl5] if one switches to the direct space using techniques developed in J52J [50] • We 
equivalently perform the renormalization procedure in the momentum space in a similar way 
of [EE]. 

The renormalization schemes corresponding to any primitively divergent contribution 
given by Table 1 have to be studied. Nevertheless, the procedure remaining similar at a 
given number of external legs, we will only treat the following cases: 

(i) iV ex t = 4 yielding ujd{G) = (first line of Table 1); 

(ii) iV cxt = 2 yielding Ud(G) = 1 (second line of Table 1). 

We emphasize that even though the above cases (and for the first situation, for a given vertex 
configuration) will be discussed, the analysis performed in the following can be carried out 
for all possible cases and affords the same conclusion. 



5.1 Renormalization of the four-point function 

Consider a four-point function subgraph G\ characterized by the first line of Table 1 and 
equipped with four external propagators. The graph is such that ggg = 0, hence has a bound- 
ary graph of the melonic type. The pattern followed by external momenta are necessarily 
of the form V4 (see Figj2|. Denoting external momenta data by p c j xt , / = 1, 2, 3, 1', 2', 3', we 
assign to each p e j Kt an external or boundary face (there are six of them). Calculations will 
be only made for the initial configuration of the V4 vertex given by ([8]). It will be obvious 
that for any other permutation, the derivations will lead to similar conclusions. 

Let Ai(G$) be the amplitude of the resulting graph. Note that the external leg indices are 
at scale ji which must be strictly smaller than i the scale index of the quasi-local subgraph 
G\- 

Using the face factorization, A±{G\) can be written (alleviating notations, henceforth p s 
refers directly to the absolute value \p s \ and the reference graph G\ will be dropped) as 

^ 4 [{^ xt }] = ^ [[]Jda e e- a t m ]{\ Y[ e -^f at)a f p T^ \ \\ e -(E*E/«*)«wU (48) 

each ae G [M~ dt , M~ de+1 ], i runs over all lines of the graph. Referring to an external line, 
we will use instead I. Hence, de can be either an external index ji or an internal one t£. In 
particular, for all / G F int in the second product, at is at scale %i 3> ji. 

Next, we single out for each independent external face the exponent otip^ , where a\ G 
[M _j! , M _j!+1 ], such that, we can write a given external face amplitude as 

e -(E« 6 /««)a/P/ xt _ e -K+"i')«/P/ xt e ~(T,ief/e^i a e)a f pf\ f^g\ 

Noting that all ae^i, in the sum, are now at scale ie 3> ji, and pj xt ~ M Ji , we perform a 
Taylor expansion for each external face as 



-(Ef 6 /^)«/P/ xt = p-(^+av) a /P/ x 



l-( V a e )a f pf t [ dte-t^ef/^^fPry (50) 
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Rewriting Q as e~^^ a{)a f p T = e -(°!+« ! ')»/pf t [i _ #ext] ; where ^ext ig the remainc ier of 
the Taylor expansion, we substitute the latter expression in the initial amplitude *4.4[{py xt }]. 
In loosely notations, the result is 



Pf i 

j "Q g-^+aj/Ja/pf j ^_ Rf + R T R f 



(51) 



/eFe X t /eF cxt /J'eFext 

Collecting the leading order contribution, one has 

e -(«!+"i')«/P/ x 



| e -(Etef a t) a fPf J j. 



/6F ir 



p/ 



n 



n »■ 

/SFext 



/6Fext 
~ext 



e -(T,eef a e) a fPf 

/eF int 



P/ ^ /eF int 



}(52) 
a /p/ 1 . 



The first factor corresponds to a vertex hooked to four external propagators since it can be 
recast in a mere form 



n da 



-aim] 



n «- 

/6Fext 



(ai+a^afpy 1 



| e_a!l ( aipl+a2P2+a3P3+m )e _a ' 2 ^ aip ' 1+a2P2+a3P3+m ^ 

7 



(53) 



e 'a 



a! 3 (aiPi+a2P2+ a 3P3+ m ) e - Q i 4 ( a iPi+ a 2P2+ a 3P3+ m ) | 



The second factor coincides with the logarithmically divergent internal contribution given 
by the power-counting theorem. Hence ^4.4[{pj xt }; 0] determines a log-divergent counterterm 
participating in the vertex renormalization. 

The remaining terms in (51) should be analyzed. Proving that the first subleading or- 
der contribution improves the power-counting will be a sufficient condition to ensure the 
convergence of the subsequent product terms. The first subleading term is of the form 



pf 



-aim] 



n «■ 

/6Fext 



(a l +a l ,)a f pj c 



(54) 



/eF cxt eef/e^i Jo /eF int 



/eFext eef/e^i 
and can be bounded by 



Pf J t 



-a e m-\ 



e -CLeef a e)afPf 

/eF int 



(55) 



where K is some constant (which could depend on the number of internal lines of the graph 
and number of external faces). Note that the integral over t yields just a factor of 0(1). Also 
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remark that, we used optimal bounds such that p^ xt < M e ^\ and, for any internal decay, 

one assumes that at < M~ l( - G i\ recalling that e(Gf) = swp t ji and i{Gf) = mi i£G ki£. The 

factor M~W G ^~ e ( G ^ guarantees that the power-counting M Wd( ^ G ^ =0 (which corresponds to 
the last sum, up to some constant) is improved and will bring the sufficient decay to ensure 
the sum of momentum attributions in the standard way of [52] • 



5.2 Renormalization of the two-point function 

We consider a graph G\ with two external legs with topology as dictated by one the cor- 
responding rows of Table 1. Note that the boundary data p^ xt , / = 1,2,3, now follow the 
pattern of a simple mass vertex. Only, the case of a two-point graph with a linear divergence 
will be treated since the subtraction of log-divergent cases can be recovered easily from the 
same analysis by restricting the Taylor expansion at less order. 

Let -4.2[{P/ xt }] be the full amplitude associated with the graph G\ with external propa- 
gators. Having defined all tools in the previous analysis of Subsection |5.1[ we perform the 
following expansion of any boundary face amplitude occurring in -4.2[{P/ xt }] a s 

l-( a l> a fPf + [-( Yl a ^ a fPTY f 1 dt(l-t)e- t{ ^f/^ ae)a f p T 

The same @ can rewritten as e'^f^^T = e -^+^') a fvT [l + Sf* +R ^ty Substituting 
this expression in the initial amplitude yields 



M{pf}} = J2 An** 

j -q e -K+^w P r][i+ y s 7 t+ E 

/S-Foxt /6-Fcxt /6-Poxt 

+ {S} xt +Rf){S}f t +Rf t ) + Yl e~^f ae)a f p f y (57) 



R f 



As expected, the leading order contribution ^42[{p} xt }; 0] is mainly of the factorized form 

ipi 



(52) with, in the present instance, the first factor appearing as 

An^ e -«n[ n 

J i /eFoxt 

= / ^n rf«i]e _a ' 1 ^ aiPl+a2P2+a3P3+m ^e _ai2 ^ aiPl+a2P2+a3P3+m ^ (58) 
J i 

and the second factor bringing a linear divergence. Clearly, this term belongs to a mass 
renor malizat ion . 

The remaining task is to prove that higher order terms improve the power-counting. Let 
us focus on the following 



v t J e 
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n «■ 

/SFext 



£ ( <»t)*fPT][ II "" ElE, °' )w ]} ( 59 ) 



/eF cxt e&f/e^i /eF iD 
which can be seen as a wave-function renormalization. Indeed, we can recompose it as 

3 



3=1 



"J^J ^ Q , ; j e - a !i( a lPi+ a 2P2+a3P3+m) e -a i2 (aipi+a 2 p2+a3P3+m) a ^ext 



P/ 



n ^ 

/6F int 



(Ef £ f a>e)a.fPf 



(60) 



where the first factor is nothing but a wave function type contribution (and renormalizes 
a s ) and the second factor is log-divergent. The latter statement is justified by the fact that 
the factors 2^e/ s /^z o.g_ brought by the derivatives are of order M~ H which annihilates the 
internal divergence degree which was linear. 

As another consequence of the decay of J2eef/e^=i ai > an y monomial in Sf of order higher 
than 2 in (57) is simply convergent. 

Next, the first order remainder term is of the form 



Pf 



-aim] 



n 

/eFcxt 
i 



-(a l +a l ,)a f pf c 



(61) 



J2 IK J2 a *) a fpTi dt(i-t)e- t{ ^w a * )a f p T 



n 

/£F int 



{J2i<zf a i) a fPf 



Using the fact that p e ^ < M e ^ G ^ whereas ag ~ M le , with i £ 3> e(Gf), for £ ^ I, we find 
the following optimal bound for the above term as 



\R 2 \ < KM-WW-'W* / []I da - 



Pf 



-aim] 



(62) 



where the integral over t yields again a factor of 0(1). We recognize, in the last sum, the 
linear divergence associated with the two-point graph at t = 0. Hence, this term is convergent 
and will ensure the summability of scale attributions. 

All the rest of terms are manifestly convergent from the fact that they involve only 
product of convergent terms. Finally, the summability of the scale attributions can be 
performed in the way introduced in 



5.3 About enlarging and reducing the space of couplings 

It is worthwhile to discuss how the (perturbative) renormalizability of the initial model can 
be extended to different models with various degrees of anisotropy between colors. Different 
classes of models are obtained by just putting a different coupling for each term or conversely 
identifying the different couplings. 
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Renormalizability is triggered by the coherent combination of three ingredients: a multi- 
scale analysis, a locality principle and a power- counting theorem. The issue raised above can 
be rephrased as "is a given class of theory with particular anisotropics stable or not under 
the RG flow?". 

Consider the renormalizable model defined by the action written in terms of Q and (J8|, 
which we call in the following discussion (I). The model has many wave function couplings, 
namely a s , and one single interaction coupling, A (the mass is always kept fixed in the next 
developments). 

In addition to the model (I), let us discuss three other basic classes of models: 

(1) The most general "anisotropic" model which has different couplings for each interac- 
tion and different wave function couplings; hence, in this model, the main parameters are 

(^e=l,2,3; a e=l,2,3); 

(2) The model with "anisotropic interactions" which has different interaction couplings 
but a single wave function coupling (X e =i,2,3, Qe=i,2,3 = a); 

(3) The "isotropic" model which has a single coupling for interactions and a single cou- 
pling for kinetic terms (A e= i ) 2,3 = K a e=i,2,3 = a )- 

Note that all these theories are continuously connected by parameter deformations. How- 
ever, the interesting question is whether these classes are stable under the renormalization 
group flow. 

All models (l)-(3) share a common feature with the model (I) studied so far: their multi- 
scale analysis and power-counting theorem are identical since those are written without the 
explicit knowledge of the coupling constants. All what is required to perform a power- 
counting is the type of propagator decay, the vertices and lines combinatorial properties and 
gluing rules and these are exactly the same for these models. One can call these theories 
"power-counting renormalizable" to this respect. 

However, as previously emphasized, the renormalizability is also about a locality principle 
that ought to be satisfied by the model. This feature can be also called "replicability" or 
"surviving property" along scales of the model. The locality principle in a just renormalizable 
model with marginal terms with n-valent vertices ensures, roughly, that any log-divergent n- 
point function with external data of the same form of a given interaction should renormalize 
the coupling of this interaction. This principle holds immediately either for a unique coupling 
or different ones for marginal terms. Indeed, one has just to decide, in any situation, to affect 
to the renormalized coupling constant of a given vertex all divergent terms with akin external 
data. Even for a unique coupling shared by several interactions, this can be explicitly done. 
Of course, from the point of view of the RG flow, the RG equations can be very involved 
(the following analysis performed at one-loop does not displays this feature but it certainly 
occurs are further loops) and each renormalized coupling is generally dependent of many (if 
not all) couplings. The key point is that they can be, at least, written in an unambiguous 
way, for one or many coupling constants for the interaction. 

The actual and only issue during the renormalization occurs potentially with the locality 
principle for relevant terms. In the present story, relevant terms are two-point functions^] 
given by the second line of Table 1 . Note that the second and third lines consist in marginal 
two-point functions and so can be simply handled by the mass renormalization. 

6 In several other contexts, they are always the most diverging ones. 
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Consider model (I) as starting point. Performing a Taylor expansion of two-point func- 



tions around their local part (in the way already introduced in Subsection 5.2) yields two 
important corrections: a 0-th order linearly divergent contribution which again contributes 
to the mass and a first order log-divergent contribution involving a s p ext . For each two-point 
function, one finds that the second kind of correction has two features: (A) it can be identi- 
fied with a term present in the Lagrangian and (B) it should contribute to the wave function 
renormalization Z s of the model. 

Note that, for theories (I) and (1) having many coupling constants a e= i ) 2,3, there is no 
ambiguity to renormalize independently each of these couplings by the above log-contribution 
of any two-point function and, so, to define Z s . This is totally independent of the presence 
of several coupling A e= i 2,3. It becomes immediate that (I) and (1) satisfy a locality principle 
for two-point functions and therefore are renormalizable. 

Let us focus on model (3) assuming that A e = A, and a e = a. Condition (A) above is 



fulfilled since (60) can be reached in any case. However, due to the fact that we have a 
single wave function coupling, Condition (B) changes drastically: It should be replaced by 
(B'): the local log-divergent contribution of the two-point function should define a unique 



wave function renormalization Z independent of the strand index s. One notes from (60), 
that the log-divergent term actually involves the strand index s (see the last line where 
Y^e.£f a /ijki a £ IS integrated). This issue can be overcome by a permutation of the vertices 
present in the theory. The latter statement is understood in the sense that given a graph 
which is not symmetric in the strands, there always exists a set of graphs (with the same 
number of vertices, lines and the same degree of divergence) identical by color permutation 
to the former. This simply holds by permutation the vertex indices V41 — > —> V^. 



Hence, if the next order correction of a two-point graph has an amplitude A' 2 (60) depending 
on s — 1, the same correction has two partners of the same form for s = 2,3. Adding 
these restores the symmetry in s and entails that the self-energy £ [p| xt , p| xt , p^ xt ] (the sum 
of all amputated 1PI two-point functions at external momenta p ext ) should be a symmetric 
function in all strands. The wave function renormalization 

Z — 1 cLcxt E [p cxXl , p c 2 xt j ^3 Xt ] \ P c f =o = l <9pS xt ^ bi xt 5 PT* > PT* ] \p c f =0 

= i--d pr i:[pT\pT\pT%T=o (63) 

should not depend on s. By expanding £ [bf^ , , b^ xt ] order by order and even though 
the latter is symmetric, the fact that Z keeps a fixed value for all derivatives d^t cannot 
be met unless all coupling coincide: A e = A. This condition clearly holds for the model (3) 
and therefore the latter is renormalizable. In contrast, the model (2) having distinct X e is 
not stable under RG flow: starting from a bare action in the class (2), one ends up with an 
effective action in the general class (1). 

Let us discuss, finally, the following peculiar partly anisotropic theory 
(4) Consider the model with two interaction couplings and two different wave function 
couplings such that Ai 7^ A2 = A3, a\ 7^ a 2 = 03; note that we could have chosen any 
symmetric situation. 

The model (4) is interesting because, we do not choose the easy case to let all a s to be 
different (which will obviously lead to renormalizability along with theories (I) and (1)), but 
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we choose instead two wave function couplings a 2 = 03 7^ a,\- This model is renormalizable 
because it lies exactly in between the two different situations so far carried out. 

Merging the sector a 2 = a 3 and following step by step the above analysis, it is direct 
to find the necessary condition A2 = A3 for having a unique wave function renormalization 
Z 2 = Z 3 . We conclude that the model (4) is renormalizable. 

The model (4) is interesting by itself for the remaining analysis on /3 function. Indeed, 
it is related to two renormalizable models which are, furthermore, asymptotically free (see 
next section): 

(4') Ai ^ \ 2 = A 3 = 0, ai ^ a 2 = a 3 ; 

(4") = Ax ± A 2 = A 3 , a 1 ^a 2 = a 3 . 

It is not excluded that there exist other renormalizable theories in this framework. The 
above analysis illustrates the richness of tensor field theories. 



6 one-loop 7-, /3- functions and RG flows 

This section starts the second part of our analysis which aims at calculating the so-called 
7- and /3-functions governing coupling constant RG flows of the model. We will restrict the 
study at one-loop since it is well-known that important properties (like asymptotic freedom 
for instance) of the model can be inferred even at this approximation. 

The U(l) 3 model described so far has a unique coupling constant associated with all 
vertices V4. It is interesting to relax that condition and study the same model in full generality 
by assigning to each vertex a different coupling constant. According to the discussion of 



Subsection 5.3, the model obtained in this way is again renormalizable. The RG equations 



for the renormalizable model (3) (from Subsection 5.3) can be inferred from this extended 
model. Meanwhile, RG equations for models (4') and (4') have to be computed with slightly 
modified method. Moreover, we introduce a symmetry factor due to the internal symmetry 
of the vertices: each coupling constants will be of the form — A e /2, e = 1, 2, 3. 

In the following, the basic quantities (the self-energy and T 4 -function) are first computed 
in full generality (i.e. within the framework of model (1)) and, then, we will particularize 
these on three reduced theories given by 

(i) [a e = a, X e = A] 

(ii) [a 2 = a 3 ^ a u A 2 , 3 = 0, Ai ^ 0] 

(iii) [a 2 = a 3 ^ a u A 2 = A 3 ^ 0, Ai = 0] 

The instance (i) corresponds to an interesting model with a unique coupling constant for 
all interactions ^ and a unique wave function coupling. It is simply the most natural 
model which can be viewed as the rank-3 analogue of the model investigated in |45j . On 
the other hand, the second and third correspond to some more drastic truncation that one 
could perform. The latter (ii) and (iii) share some but not all features of the former hence 
are different. 
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6.1 7" an d /3- functions 



In full generality (working in the model (1) Subsection 5.3), considering three wave function 
renormalizations ^=1,2,3, each of one with respect to each of the strands, the field strength 
can be modified as follows: 

<p —> (Z&Za)* <p, (64) 
so that, after renormalization, the wave function couplings satisfy the equations 

aT = a e (J^j 3 , 6 = 1,2,3, e ^ e, (65) 

For a rank three tensor model as the one we are presently dealing with, we define 

1 



AS 



(66) 
3=0 



where the self-energy £(&i, b 2 , b 3 ) is the sum of the amputated one particle irreducible (1PI) 
amplitudes of the two-point correlation function truncated at one-loop that we denote 

&2, 63) = ( </%fc 2 &3 ^6x6263 )ipi- ( 67 ) 

Our initial goal is to compute at one-loop the dynamics of the effective couplings a e governed 
by the 7 e -functions encoded in (65). In a second stage, we will compute, still at one-loop, 
the mass /3 m -function given by the expression 

m „„ = m-EtO.O.IJ) m 
(Z 1 Z 2 Z 3 ,) i 



I 3 



In the final step, we will study the dynamics of constant couplings A e , governed by the 
/3 e -functions ciphered by the following equations 

r = _WW) ! f = 123 m 

(ZiZ 2 Z 3 ) 3 

where r4, e (&i, b 2 , 63, &i, b' 2 , b' 3 ) is the amputated 1PI four-point function. 



Let us recall that the amputated 1PI four-point functions occurring in (69) read 

r4,l(&l, b 2 , 63, b[, b' 2 , 63) = ( <£> 6l & 2 & 3 Vb'^bs W&Vs ^6i& 2 & 3 )h 



IP/' 



P^&l, b 2 , 63, b[, b' 2 , 63) = ( ^6^263 <£>bi& 2 & 3 ^6^263 ^6263 )lP/> 

r^C&l, & 2, &3, b' 2 , 63) = ( </?6i6 2 & 3 ^bi6 2 6 3 ^6ib' 2 6 3 ^6^,63 )lP/> ( 70 ) 

where external indices, even though repeated, are not summed and follow the pattern of the 
vertices of the model. This is justified by the renormalization prescription. 

Considering the reduced cases (ii) and (iii), there are potentially two wave function 
renormalizations Z e= i >2 . Nevertheless, the cancellation of one or many couplings may have 
drastic consequences on the way that the /3-function equations have to be written. We will 
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deal with these after computing the self-energy as well as the r 4 -function in full generality 
and then putting to zero some of the contributions. 

For €,£,<= = 1, 2, 3, and pairwise distinct, we introduce the following formal sums: 



S e (b) 


:= ^ l/(a e \b\+ 0,^1+ ai\p 2 \+m), 

Pl,P2 6S 




(71) 


S' e (b,b>) 


:= '£l/{a e \b\+a i \lf\ + a e \p\+™)> 
pes 




(72) 


S e 


:= ^2 1 /( a e\pi\+ az\p 2 \+m) 2 , 
Pi,p 2 es 




(73) 


s[ 


= ^l/(a 2 (| Pl | + |p 2 |)+m) 2 , S> 2 = 

Pl,P2 


^ l/(ai|pi + a 2 \p 2 \ + m) 2 , 

Pl,P2 


(74) 




:= ^ 1 /[( a e\ b \ + a e\Pi\ + a'M + m)(a t 


\b'\ + ai\px\ + ai\p 2 \ + m)) 


(75) 



PI -P2 



Note that S t is linearly divergent whereas S' e , S e , S' 12 and S t are logarithmically divergent. 
Another important fact to notice is S e (0, 0) = S t . 
Let us prove the following proposition 

Lemma 5. At one-loop, the self-energy and wave function renormalizations are given by 
S(6i, b 2 , b 3 ) = X e S 6 (b 6 ) - XtSHh, b 3 ) - X 2 S 2 (h, b 3 ) - X 3 S' 3 (h, b 2 ) + 0(A 2 ), 

c 

Z e = l-X t S t + 0(X 2 ), 6 = 1,2,3, (76) 
where 0(X 2 ) is a O-function of any quadratic products of any coupling constants Ai^- 
Proof. We can first evaluate the self-energy mj\ as 



E(h,b 2 ,b3) =J2 K 3 A 5(Kb 2 ,h), (77) 
Q 

where Q runs over two-point 1PI graphs with amplitude Ag(pi, b 2 ,b 3 ) and corresponding 
combinatorial weight Kg. The latter is known to be the number of Wick contractions given 
rise to Q. At one-loop, only the tadpole graphs T 1; T 2 , T 3 (such that J2j9j = 0) ^2 
and T3 (such that J2j9j = 1) (Fig|9] provides T\ and T[, the remaining configurations can 
be easily recovered by color permutations corresponding to the other configurations of the 
vertex ip ) contribute to (77) with the combinatorial factors 

K Tt = 2 = K n , e = l,2,3. (78) 

One gets the amplitude for each tadpole T e 

A Te (b e ) = -^S e (b e ), (79) 



where S e (b e ) is given by (71). Meanwhile, for T £ ', we have the amplitude 

A^(Ml) = (80) 
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Figure 9: Tadpoles for the first type of vertex. 



with S' e (be,bt) given as (72). At first order, we obtain the self-energy as the sum of these 
contributions as 

£(&!, h, h) = ~J2 X * S M ~ Al5 i( &2 > 6 ») ~ ^(Ms) - X 3 S' 3 (b h 6a) + 0(A 2 ). (81) 

It should be emphasized that not all contributions of 2(61,62,63) have to be taken into 
account for the wave function renormalizations (66). Only those leading to a log-divergent 
behavior after differentiation have to be considered. A quick inspection shows that S' e (b,b') 
are log-divergent hence should be only considered for the mass renormalization. The wave 
function renormalizations (66) can be finally expressed as 

Z e = l-X e S e + 0(X 2 ). 



(82) 



where S e is given by ( 73 ) . □ 
We are in position to compute the dynamics of the renormalized wave coupling constants 
a v e en and of the renormalized mass m ren . In the same previous notations, the following 
statement holds: 

Theorem 4. At first order, the renormalized wave function couplings and mass satisfy, 
respectively, 



a™ = a f 



m ren _ m 



3 

+ ^A e S e (0) + O(A 2 ,AlnA). 

6=1 



+ 0(A 2 ), e = l,2,3, e^e, e^I^e,(83) 



for momentum cut-off A. 



Proof. Concerning the first statement (83), by simple index permutations, all cases can be 
easily deduced from e = 1 hence we will focus only on this situation. Using Lemma [5] and 
remaining at first order in the constant couplings A e , we have: 



Gti 



1 -2A 1 g 1 + Q(A 2 ) 
1 - (X 2 S 2 + A3S3) + 0(A 2 



cii 



+ 0(A 2 



(85) 



Focusing on the renormalized mass equation (68), Lemma |5| allows us to write 

m + E £ KS e (0) + X 1 S[{0, 0) + A 2 S 2 (0, 0) + A 3 5 3 (0, 0) + 0(A 2 ) 



m 



1-UX 1 S 1 + X 2 S 2 + X 3 S 3 ) + 0(X 2 ) 
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= m + J2 KSt{0) + XiS[{0, 0) + X 2 S' 2 (0, 0) + A 3 ^(0, 0) 

e 

+ 1 (XiSx + X 2 S 2 + A 3 S 3 ) + 0(A 2 ). (86) 

Neglecting the subleading divergences compared to the linear divergence of «S e (0), one is led 
to (84) after having introduced a cut-off A in the sums S' e (0, 0) and S e . □ 
Discussion. The 7-functions of the model can be defined by restricting the space of the 
couplings parameters to a smaller subspace. In fact, the reduction will be performed on the 
space of six couplings Ai, 2 ,3 and 01,2,3- 

Many cases may occur by collapsing couplings. The first case is straightforward 



(i) We set a e = a in (83). This enforces Ai = A2 = A3 and induces a unique equation: 

a ren = a. (87) 

The latter means that the wave function coupling becomes stationary in the UV. In 
other words, the 7-function for this case is trivial: 

7 = 0. 



Meanwhile, from (84), one infers 

Pm = -3. (89) 

(ii) Otherwise, setting A2 = A3 = but Ai 7^ yields necessarily that a\ 7^ a 2 . We need to 
re-evaluate the self-energy from (81) and we obtain the only non trivial wave function 
renormalization 

ZP = i-A 1 S , i + 0(A?). (90) 



where S[ is given by (74). The equations of the couplings find some modified form: 

aS' )rcn = ai, 4' )rcn = a 2 ±- = a 2 (l + A^) + 0(A?) . (91) 
The 7-functions at one-loop in this restricted space are given by 

1? = 0, lP = -1. (92) 

Assuming that the sign of Ai is positive, a\ has a stable value in the UV, whereas a 2 
flows towards a vanishing value. Under the same assumptions, we get for the mass, 
the renormalized mass equation and f3 m function 

m ( ' )ren =m + Ai5i(0) + O(A?,AilnA), P$ = -1, (93) 

where <Si(0) can be inferred from <Si(0) after identifying a 2 = 0,3. 

(iii) Setting A 2 = A 3 7^ 0, Ai = and a x 7^ a 2 and computing the self-energy, one has the 
unique non trivial wave function renormalization 

ZP = l-X 2 S' 2 + 0(X 2 2 ), (94) 
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where S' 2 is provided by (74). In this restricted space, the coupling equations are given 
by 



(") ren 



02, 



(")ren 
1 



1 



zp 



a 1 (l + X 2 S' 2 ) + 0(X 2 2 ) 



(95) 



,(") 



■1. 



yielding the 7-functions at one-loop 

iP = 0, 7i 
The same assumptions give for the mass 

Thus cases (ii) and (hi) share similar properties. 

Theorem 5. The renormalizable coupling constants associated with X e , for e = 1, 2, 3, satisfy 

2 



(96) 



(97) 



a; 



K + K 



—XfSf 



AiSi + A25*2 + A3S3 



0(A d 



(9S 



where 0(A 3 ) zs a O- function of any cubic product of all couplings Aj.,2,3- 
Proof. The one-loop 1PI four-point functions can be written as 

r 4 ,e(oi, 62, 63, b[, b' 2 , 63) = K Q AK&i, 62, 63, & 2 > & s) 



(99) 



where ^ is a four-point 1PI graph with the topology as required by the renormalization, with 
amplitude Ag(bi, b 2 , 63, b[, b' 2 , b' 3 ) and combinatorial weight Kg. 

Note that, interestingly, the "mixed" graphs obtained by gluing vertices of different kind 
do not contribute to the effective coupling constants. Indeed, although these graphs could 
form melons, they all possess two lines and one face giving a convergent power-counting 
M -1 . Another way to agree with this fact is to observe that the boundary graph (even 
though melonic) of any of them is disconnected with Cqq = 2. As a consequence, in the 
following developments, cross terms involving product of coupling constants A e A e /, e 7^ e', are 
inexistent. Hence, a unique graph F € contribute to r 4 e (Fig, 10 displays such a contribution 
for e = 1, the remaining can be deduced by permutations) and its combinatorial factor is 
always 

K Fe = 2-2 -2. (100) 
Given an e, the amplitude of each graph is such that 



2 2 2 



(101) 



where S e (b e ,b' e ) is given by (75). 

The amputated 1PI four-point functions 1^=1^,3 defined by (70) and (99) can be evalu- 
ated at one-loop as 



M&i, & 2 , 63, &;, b' 2 , b' 3 ) = -x e + x 2 e s t (b € , K) + o(x 3 



(102) 
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Figure 10: The unique melonic 1-loop four-point function with 
melonic connected boundary for the first type of vertex. 



where 0(A 3 ) is a O-function of any cubic product of coupling constants. Therefore, at low 
external momenta, they reduce to 

r 4)6 (0, 0, 0, 0, 0, 0) = -A e + X 2 € S € + 0(A 3 ). (103) 



Having all required ingredients, we are in position to evaluate the ratios (69). Using Lemma 



5jand (103), by direct algebra, one comes to 

-A e + X 2 e S e + 0(A 3 ) 



{Yl e (i-x e s € + o(x*))} 



Ae + X e (^ - AeSe + \ (Ai^I + X 2 S 2 + A 3 S 3 ) j + 0(A 3 ) 



(104) 

□ 



which achieves the proof. 

Discussion. We can once again discuss the merging of coupling constants Ai 2,3 into specific 
varieties in order to deduce the /3-functions. 



(i) We merge all couplings such that a £ = a and X e = A, such that S € (73) becomes 

S= VKbi| + b 2 |)+m)] 2 , 



(105) 



then we have from (98) 



X mn = X + X 2 S + 0{X 3 ). (106) 

such that the /3-function of the model with single wave function renormalization and 
single coupling constant is given by 

P = -l. (107) 
This model is therefore asymptotically free. 

(ii) Assuming that a\ 7^ a 2 = a 3 , and Ai 7^ = A2,3, there is a unique equation for the 
coupling constant Ai expressed as 



(') ren 



r C) 

4,1 



From this, we infer 



-A! + A^ + Q(A 3 ))) 
(1-A 1 S 1 + 0(A?)) 2 



Pi 



(') 



Ai + A^i + 0(A;). (108) 



(109) 



Consequently, the model is also asymptotically free. 
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(iii) Last, we set A 2 = A 3 7^ 0, Ai = and ai 7^ a 2 . The coupling constant equation can be 
written 



A 



r (") 

(") rcn _ 1 4,2 



{z py- 



(110) 



and yields in the same way as previously done in (108) the one-loop /3-function 



111) 



showing that this theory is asymptotically free. 



6.2 RG flows 



We restrict the RG flow equations to the particular situations discussed so far. It turns out 
that for these cases, explicit solutions are affordable at this truncation. Cases (i) and (ii) 
will be only discussed for simplicity. 

Case (i): A e = A and a t = a. This case yields the coupling equations (87) and (106) of the 
form: 



a, 



A r 



A + A 2 S + 0(A 3 



(112) 



Hence, the wave function renormalization is stable whereas the coupling constant satisfies 
the discrete RG equation (truncated at first order) 



Aj-i — Aj + A 2 Si, 



113) 



where Si is given from S ( 105 ) by restricting the sum to high momenta \pi | + |p 2 | G [M l 1 , M l ] 
such that (see Appendix Bjfor more details) 

S i =^K M + 0(M- i ), (114) 



with km a log-divergent term in M which does not depend on the scale i. We get the 
solution, using = a, 



d\i 



km di, 



a 2 A, 



4«: M (A - i)K 



115) 



where we have introduced the coupling value X uv at the cut-off scale A. The RG flow of A is 



pictured in Fig 11 

Case (ii): A 2 = A3 = 7^ Ai and a 2 = 03 7^ a\. The RG flow equations can be obtained 
from ([9il and Jl08|) as 



a 
A? 



rcn 
1 

rcn 



ai, 



A! + A 2 ^ + 0(A?; 



a 2 (l + A!^) + 0(A 2 ), 



;ii6) 



with S[ given by (74). The same RG flow equations (truncated at one loop corrections) 
translate in discrete version at scale i as 



01; i-1 — Q>1; ii a 2; i-1 — «2; i [l + Ai; j S*^ j] , Ai ; j_i — Ai ; j + A 2 , j S*^. 



117) 
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a 2 = 1/2, 1 ; k = Log 2 ; A uv = 1CT 4 ; A = 1000 



X 




200 400 600 800 1000 

Figure 11: RG flow of the coupling constant A running to a trivial fixed point 
for two values a = 1/2 (top) and a = 1 (bottom) and for a UV cut-off A = 1000 
and other parameters fixed such that X uv = 10 and n = logM = log 2. 



Summing S[ only on high momenta \pi\ + \pz\ G [M ,_1 ,M'], we have = 4/o|. ^ka/ + 
0(M - *). Focusing of the second and third expressions, these discrete equations can be put, 
at first order, in an equivalent differential form 



da2- ■ 
di 

d\i- 
di 



a 2; 
1 



a 



\ 2 

A x . i Km 



which induces that 



2: i 



for some constant K. We substitute K in (118) and get 

— AkmK 2 di 



dX\. 



Ai(i) = A^ + 4k m ^ (A-i). 



(118) 



119) 



(120) 



Thus, one can readily obtain a%{i) by combining (120) and (119). We see that the flow of 
this model actually differs from the previous one. Both couplings Ai and 02 go linearly to 
their UV fixed values. 



7 Conclusion 

We have shown that a rank 3 tensor model on U(l) 3 , analogue of the previous model identi- 
fied in , is perturbatively renormalizable at all orders of perturbation theory. The proof of 
this statement relies on both combinatorics and the colored model properties serving as un- 
derpinning of the present class of tensor models. The interplay of three ingredients, namely 
a slice decomposition, a power-counting theorem and a generalized locality principle has 
guided us towards a proof of renormalizability of this model. We have also introduced differ- 
ent wave function couplings for each strand hence giving them a different dynamics as well as 
different interaction coupling constants. This definitely enlarges the class of renormalizable 
models of the kind. As it should be also emphasized, the renormalization has been performed 
in the momentum space and not in the direct space. This particular basis has allowed us 
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to identify an analogue theorem for the momentum routine useful for renormalization in 
ordinary quantum field theory [52J . The second part of this contribution was devoted to the 
computations of the 7-, /3-functions and the corresponding RG flow of couplings constants 
in the theory. There exist two classes of underlying models turn out to be asymptotically 
free in the UV limit: (a) The model obtained by merging all coupling constants, such that 
^1,2,3 = A and ai )2 ,3 = a; (b) the class of models defined by A 2 3 = 7^ Ai and a,\ 7^ a 2 = a 3 
and the symmetric obtained by color permutation. 

The model introduced here and its analogue developed in jl5] claimed to be simplified 
but possible models for quantum gravity using the "tensor device" |46j . What we have shown 
here, by exhibiting such a renormalizable and asymptotically free (in the UV) model, is that 
the tensor approach for quantum gravity is a promising line of research for those who believe 
that gravity should be described by a renormalizable quantum field theory. 

Future prospects are can be performed in different directions. Nonlocal tensor models 
over copies of M can be treated provided one introduces an IR regulator. For instance, 
the Mehler propagator l/(p 2 + x 2 ) (on each strand) could be an interesting proposal. The 
resulting models should be called tensor extension of the Grosse-Wulkenhaar model [IE] 
in noncommutative field theory. Another important question is whether or not such an 
analysis could be handled for GFT models on copies of £77(2). A priori, since the tensor 
1/N expansion is also valid for GFT models, the answer is yes. Nevertheless, the level of 
difficulty is much higher is that situation. The renormalization will definitely involve both 
£77(2) recoupling theory and Taylor expansions of graphs around their local parts. These 
later entail generalized saddle point analysis which are far to be easy for general graphs at 
any order of perturbation theory [31] . Besides, going deeper in the analysis of the class of 
models presented here, an important point would be to investigate if the rank four model 
shares this important property of asymptotic freedom. Finally, interesting though more 
abstract questions can be also addressed at this stage. For instance, the determination of 
topological polynomials associated with the graphs of these models using an extended (i.e. 
tensor) parametric representation of graph amplitudes in the spirit of [5H] could be a fruitful 
line of investigation. 

Appendix 



A Proof of Lemma 

A.l Proof of Lemma [3] 

This appendix provides the proof of Lemma [3] using a sequence of Ofc-dipole contractions. 
First of all, we need to define a central notion which is the following: 

Definition 3 (Pairs). Let Q co i or be a colored graph and consider a Ok-dipole inside Q co i or . A 
"pair" is a couple of colors (ab), a,b = 0,1,2,3. A pair is called "outer" if the two colors 
are external to the dipole. A pair which has one color inside the dipole and one out is called 
a "mixed" pair. A pair with two colors inside the dipole is called an "inner" pair. A pair 
(ab) is said to belong to a jacket if the pair is one of the four adjacent (consecutive) pairs in 
the jacket cycle (Oxyz) . 



3 and divergent graph classification 
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Figure 12: Different types of pairs highlighted for each Ofc-dipoles: 
External pairs (in red), internal pairs (in green) and mixed pairs (in 
black) . 



The total number of pairs is always 6 and the number of mixed pairs is at least 3. Fig, 12 
indicates the different kinds of pairs. 

Contracting a dipole leads to different cases and we urge to further classify the different 
pairs in a more concise manner. An outer pair is said of type A, or disconnected by the 
dipole contraction if the half-strands at each corner on the left and on the right of the dipole 
belong to two different connected components of the graph after the contraction. Otherwise, 
we call it a "special" pair. A special pair can be single-faced if the two corners belong to 
the same face of the graph, or double-faced if the two corners belong to two different faces 
of the graph. Any type A outer pair should be single-faced at the beginning. We therefore 
classify outer pairs in the following three types: 

• Type A outer pairs are single-faced, 

• Type B outer pairs are single-faced, 

• Type C outer pairs are double-faced. 

Transverse pairs are mixed pairs which do not change their number of faces under con- 
traction. Inner pairs have one face less after contraction. Type A and B outer pairs have 
one face more after contraction whereas type C external pairs have one face less after the 
contraction. In summary, for any jacket J, we can relate the number of faces before and 
after contraction by the formula 

~ F J = \ A j\ + \ B 7\ ~ ~ \h\ , (A.l) 

where \Xj\ is the number of pairs of type X=A,B,C, in the jacket, and is the number of 
inner faces. 

Lemma [3] states that, for a single jacket J, we have 

gj -g~>0. (A.2) 

We emphasize that the proof of this statement which will be given here slightly differs from 
the proof of a related statement as found in [45J. Indeed, in the latter reference, the sum 



over jackets is performed at the level of (A.2) and leads directly to a bound on ^Agj — gj,) 
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ABC 
Figure 13: 00-dipoles configurations. 

However, as was mentioned [15] and was computed in a earlier version of the same paper, 



the relation (A. 2) is a stronger statement than the result on the sum of jackets and therefore 



worth to be communicated. 

Proof of Lemma [3} We first translate the difference between genera of jackets before and 
after contraction in terms of the basic combinatorial elements: 



2 - 2gj 
2c - 2gj, 

(9J-9J>) 



V 

V 
1 



L + Fj, 
-L' + Fj, 



(V-2)-(L-4)+Fj, 



[{Fj, - Fj) + 2 - 2(c' - 1)] 



(A.3) 



The sign of gj — gj, reduces to that of (Fj, — Fj) + 2 — 2(d — 1). The latter quantity proves to 
be always positive jacket by jacket. To prove this claim, we perform the following systematic 
analysis on the types of 0/c-dipole contractions. 

• lrst Case: 00-dipole contraction. This case is defined by a contraction of a unique 
internal line with color 0. In any initial colored graph, there are three mixed pairs, three 
outer pairs and no inner pair. Meanwhile, each jacket contains two mixed and two outer 
pairs. 

- lrst subcase d = 3 (Fig, ISA). This situation happens if the initial graph has a connected 



two-point subgraph on each line 1, 2, 3. Consequently, all three outer pairs must be of type 
A. Hence, the 00-dipole contraction yields for all jackets: 



(F F) + 2-2(c'-l) = 0. 



(A.4) 



2 (FiglSB). This case happens if we have one connected two-point 



so that (46 ) is true. 

- 2nd subcase d = 

functions plus one connected four-point function on four half-lines hooked to the dipole. In 
that case, we have two corner pairs of type A and one special pair of type B or C. It can be 
seen that 

!j) + 2-2(c'-1) = 2-2|Cj|. (A.5) 

Hence in all cases, ( f46~| ) 

- 3rd subcase d 



(Fj, 
is true. 



1 (Fig.ISC). Contracting the dipole gives d — c = 0. This can happen 



when a connected six-point function is hooked to the dipole. There no type A corners and 
three special pairs. Each jackets contains at most two special pairs. Contracting the dipole 
yields 



Fj) + 2 



2(d 



2|Cj|, 



(A.6) 



such that (46) is true. 
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A 



B 



Figure 14: A jacket (0123) which contains the pair (03) defining the 
dipole (A) and another jacket (0123) from another dipole which does 
not contain the same pair. 



• 2nd Case: 01-dipole contraction. There is one inner, four mixed and one outer pairs. 
A distinction should be made between inner- adjacent jackets, for which the inner pair of 



the dipole belongs to the jacket (see Fig 14 k, for instance, the pair (03) belongs both to 



the jacket (0123) and the dipole) from non-inner- adjacent jackets for which it does not (see 
Fig, 14B, for instance, the pair (02) is not in the jacket (0123) but defines the dipole). We 



count two subcases: 

- lrst subcase d = 2. In that case, the unique outer pair is of type A. An inner- adjacent 
jacket has two mixed pairs, one inner pair and one pair of type A. A non-inner- adjacent 
jacket has four mixed pairs, no inner pair and looses the type A pair. One finds, in all cases, 



(Fj> 



Fj) + 2 



2(c'-l) =4-4 = 0, 



(A.7) 



and (46) is verified. 

- 2nd subcase d = 1. This case reduces to the presence of one special pair. An inner- adjacent 
jacket has two mixed pairs, one inner pair and one special. A non- inner-adjacent jacket has 
four mixed pairs, no inner pair and no special pair. In all cases, 



(Fj, - Fj) + 3 - 2{d - 1) = 2 - 2\Cj\ 



(A.8) 



and (46) is true. 



• 3rd Case: 02-dipole contraction. There are three inner pairs, three mixed pairs and 
no outer pair. Only the case d = 1 occurs and one gets 



{Fj, - Fj) + 2 - 2(c' - 1) 



0. 



(A.9) 



hence, (46) is satisfied. 



• 4th Case: 03-dipole contraction. This case is straightforward and corresponds to the 
contraction of a vacuum connected component (with two vertices and four lines). There are 
six inner pairs, each jacket containing four of these pairs and we have (d — 1) = — 1. Thus 



(Fj - Fj) + 2 - 2(c' - 1) = . 



Hence (46) is again valid. 



(A.10) 



□ 
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A. 2 Classification of divergent graphs 

We extract from the divergence degree data all graphs which could lead to uj<i(G) > and 



aim at proving that Table 1 of Subsection |4.3| exhausts all possibilities. 

Let us define the integer P(G) = (Cag — 1) + V 2 + | [^cxt _ 4]. We reformulate Lemma 
[4] in the following terms: 

Ud{G) = -[Y J 9j-9dg]-P{G) < -2g dg -P(G), (A.ll) 
j 

g dg = and ^ gj > u d (G) < -1 - P(G) ■ (A.12) 

j 

We now seek the list of graphs with Ud(G) > which are those which should be renormalized. 
Case iV ext > 4: Considering N ext > 6 leads to P(G) > 1- Thus, u(G) < —1 and the graph 
amplitude is simply converging. 

Case N cxt = 4: Under this condition, one has P(G) = (Cog — 1) + ^2- The divergence degree 
is at most zero. The significant case occurs when oj(G) =0 and for that one must have 

Cgg = l, g 9 g = Y j gj a =^ V 2 = 0. (A.13) 

Ja 

Case A'cxt = 2: Since Cag < N ext /2, we should have Cag = 1 and, therefore, P(G) = V 2 — 1. 
Furthermore, the only possible boundary colored graph made with two vertices is unique 
and does not have any genus gag = 0. We conclude that the divergence degree is at most 1. 
One may have u(G) = 1 if P(G) = — 1, and only if 

C dg = l, g d g = J29j = 0, V 2 = 0. (A.14) 

j 

Besides, u>(G) = can happen in two cases: (1) P(G) = 0, in which case, we have 

Cag = l, g dg = J29j = Q, V 2 = l (A.15) 
j 

or (2) P(G) = -1 and, so, 

Cgg = l, g d g = 0, Y.9J =1 > V * = " ( A ' 16 ) 



In summary, the divergent graphs are determined by the Table 1 of Subsection 4.3 



B Formal sum approximations 

We start by proving the claim ( |114[ ) by giving the relevant part of the formal sum 
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' (P+1) 1 



s t = y ^ - ; l = — ^— + (m- +i + a + 1) - ^ (0) (m' + a + 1) 

-(A - 1) (ikP + A) - (M i+1 + A + 1)) (B.17) 



where A is a constant (which is set ^ in the text) and ip^ n \z) is the so-called polygamma func- 
tion or n th derivative of lnT[z]. There are well-known relations satisfied by the polygamma 
functions and are given by 

+ z ) = ^ m \z) + (-l) m m! z- {m+1 \ 
tp(°\z) ~ In z — - — h . . . , z — > oo, |argz|<7r, 

(n — 1)! n! 

+ 



2 — > oo, I argz| < 7r, n > 1,(B.18) 



so that the divergent part of (B.17) as M l ^> 1 can be given as (neglecting all factors of the 



form constant of convergent factors of the form 1/M % or constants) 

Si ~ K[(i + 1) In M -ikiM] = k\uM (B.19) 
hence this does not depend on the scale i. 
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